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SUMMARY 


A  numerical  analysis  technique  has  been  developed  to  compute  the 
fields  inside  a  three-dimensional  arbitrarily-shaped  heterogeneous 
dielectric  or  biological  body  inside  a  rectangular  waveguide.  The 
numerical  computation  involves  a  dyadic  Green's  function  containing  a 
double  infinite  series,  which  is  evaluated  by  a  partial  summation 
technique.  The  numerical  results  compared  very  well  with  reflection  and 
thermographic  measurements,  which  were  conducted  at  Walter  Reed  Army  j 

Institute  of  Research.  This  analytical  approach  should  be  an  efficient  j 

tool  in  guiding  the  design  of  microwave  applicators  used  for  rapid  ■ 

j 

inactivation  of  enzymes.  In  addition,  the  accomplishment  represents  a 

significant  progress  in  microwave  and  electromagnetic  theory.  j 
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SECTION  I 
INTRODUCTION 


The  role  of  cyclic  adenosine  monophosphate  (cyclic  AMP)  as  an 
intracellular  mediator  of  the  action  of  a  number  of  hormones  has  been 
well  established  [1-3].  In  addition,  cyclic  AMP  is  also  considered 
important  in  the  function  of  the  central  nervous  system  [4,3].  The 
distribution  of  cyclic  AMP  in  the  brain  can  therefore  be  used  as  an 
important  tool  in  neurochemical  research  such  as  in  the  evaluation  of 
the  effects  of  hormones  or  drugs  in  the  central  nervous  system.  In 
order  to  determine  the  distribution  of  cyclic  AMP,  however,  it  is 
necessary  to  inactivate  rapidly  two  enzymes,  adenylate  cyclase  (AC)  and 
phosphodiesterase  (POE),  which  produce  and  degrade  cyclic  AMP,  if  left 
active.  Otherwise,  the  level  of  cyclic  AMP  concentration  to  be  studied 
may  be  distorted  in  a  few  seconds  after  the  animal  is  sacrificed  [6-8]. 

There  are  two  types  of  widely-used  inactivation  techniques  — 
liquid  nitrogen  freezing  and  microwave  heating.  The  disadvantages  of 
liquid  nitrogen  freezing  include  the  slowness  of  the  process,  the 
nonuniformity  of  the  freezing  pattern  throughout  the  brain  [9,10],  and 
the  inconvenience  in  post-mortem  dissection  at  freezing  temperatures. 
Microwave  heating  is  a  promising  approach  based  on  the  principle  that 
cyclic  AMP  is  a  relatively  heat-stable  substance  [11]  while  the  AC  and 
PDF,  enzymes  are  heat  labile  and  denature  irreversibly  at  temperatures  in 
the  range  of  65°-90°  C.  With  microwave  heating,  the  inactivation  of 
enyzmes  can  often  be  achieved  with  exposure  times  on  the  order  of  1-2 
seconds  or  less,  permitting  subsequent  required  dissection  of  the  brain 
at  room  temperature. 

Microwave  inactivation  can  be  carried  out  by  an  open  or  closed 
system.  In  the  open  system,  a  plane  wave  is  employed  to  illuminate  the 
animal  without  the  direct  interference  of  enclosing  conducting 
boundaries.  The  disadvantages  of  the  open  system  include  a  large  power 
source  required  and  the  necessity  of  containment  of  radiated  power  with 
a  conducting  screen  or  an  anechoic  chamber.  The  closed  system,  such  as 


a  rectangular  waveguide,  is  highly  efficient  in  power  usage  and  has 
negligible  leakage  power  which  is  desirable  as  far  as  radiation  hazard 
and  man-made  noise  are  concerned.  As  a  result,  recent  research  in 
microwave  inactivation  techniques  has  been  concentrated  in  the  closed 
system  type,  such  as  the  rectangular  waveguide  applicator  [12].  Lenox, 
et  al.,  [12]  have  shown  that  using  a  waveguide  applicator,  modified  from 
a  WR430  waveguide,  only  2.8  seconds  are  needed  to  inactivate  the  brain 

enzymes  in  a  325-gram  rat  at  2450  MHz  with  a  source  of  3500  watts. 

2 

Microwave  leakage  is  less  than  5  mW /cm  at  a  distance  of  10  cm  from  the 
applicator  surface. 

Microwave  heating  can  be  considered  as  being  generated  by  a 
continuously  distributed  equivalent  source  throughout  the  whole 
biological  body,  and  is  therefore  potentially  more  uniform  than  the 
usual  exterior  heating  produced  by  localized  sources.  The  uniformity  of 
temperature  in  a  biological  body  under  microwave  excitation  is 
determined  by  the  degree  of  uniformity  of  the  equivalent  source 
distribution,  its  thermal  conductivity ,  and  the  time  of  exposure.  For 
the  case  of  rapid  enzyme  inactivation,  the  heating  pattern  is  roughly 
proportional  to  the  equivalent  source  distribution.  Therefore, 
nonuniform  equivalent  source  distribution  generates  "hot  spots"  in  the 
brain  tissue,  which  may  distort  the  AMP  distribution  to  be  studied. 
Difficulties  in  obtaining  uniform  heating  patterns  have  long  been 
recognized  at  Walter  Reed  Army  Institute  of  Research  (WRAIR)  in  its 
continued  research  effort  in  microwave  enzyme  inactivation  [12]. 

The  physical  problem  of  waveguide  enzyme  inactivation  is 
considerably  more  complex  than  the  open  inactivation  system  employing  a 
plane  wave.  The  physical  parameters  involved  in  the  process  are 
numerous,  making  it  extremely  difficult  and  expensive  to  design  a 
waveguide  inactivation  system  by  a  purely  experimental  or  trial-and- 
error  approach.  Under  these  circumstances,  the  idea  of  guiding  the 
design  of  the  microwave  applicator  by  an  analytical  approach  appears 
attractive,  and  the  present  research  grant  was  awarded  to  Georgia  Tech 
to  investigate  the  potentiality  of  this  approach. 

The  theoretical  approach  taken  at  Georgia  Tech  was  to  develop  a 
computer  algorithm  to  compute  the  distribution  of  dissipated  power 
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inside  a  three-dimensional  arbitrarily-shaped  heterogeneous  dielectric, 
body  inside  a  rectangualr  waveguide  with  matched  termination.  The 
theoretical  work  supported  by  this  grant  has  been  partially  documented 
in  a  paper  [13]  which  discussed  the  theory  and  its  numerical  validation. 
This  analytical  technique  was  then  checked  against  thermographic  data 
generated  at  WRAIR.  Satisfactory  agreements  between  the  theoretical 
results  and  the  thermographic  data  were  observed. 

The  results  of  this  research  have  demonstrated  that  the  analytical 
approach  taken  at  Georgia  Tech  is  a  powerful  tool  which  should  help 
significantly  future  development  of  a  waveguide  enzyme  inactivation 
system  to  achieve  rapid  and  uniform  heating  in  the  brain  tissue.  The 
technical  accomplishments  are  discussed  in  the  following  sections 
together  with  recommendations  for  further  research  in  this  area. 


SECTION  II 


THEORY 


The  problem  Co  be  considered  is  shown  in  Figure  1,  in  which  a  three- 
dimensional,  arbitrarily-shaped  heterogeneous  biological  body  is 
electromagnetically  illuminated  in  a  rectangular  waveguide.  It  is 
desired  to  determine  the  temperature  distribution  in  the  biological  body 
as  a  function  of  space  and  time.  The  analysis  involves  two  steps, 
namely,  the  computation  of  the  distribution  of  power  dissipation  and  the 
calculation  of  the  resulting  temperature  distributions.  The  first  step 
involves  the  solution  of  an  electromagnetic  boundary-value  problem.  The 
second  step  involves  the  solution  of  the  heat  conduction  problem  with 
the  heat  source  represented  by  an  equivalent  distributed  source  of 
dissipated  power  density. 

Although  numerous  publications  have  appeared  in  the  literature  on 
the  analysis  of  waveguide  discontinuity  problems,  not  a  single  truly 
three-dimensional  solution  has  been  published,  according  to  a  review 
paper  by  Silvester  in  1974  [14],  In  the  related  cavity  resonator 
problem,  finite-element  [15]  and  TLM  [16]  methods  have  been  employed  to 
compute  the  resonant  frequency  and  field  distribution  when  dielectric 
blocks  were  placed  inside  the  cavities.  Thus,  the  theoretical  work 
performed  in  the  present  grant  represents  a  significant  advance  in 
microwave  and  electromagnetic  theory.  Most  of  the  details  of  the 
numerical  analysis  have  been  discussed  in  a  publication  [13]  which  is 
included  in  this  report  as  an  appendix.  However,  the  treatment  of  the 
singular  problem  in  the  dyadic  Green's  function  was  only  briefly 
mentioned  in  Reference  13.  Because  of  the  controversy  in  the  literature 
on  the  singularity  problem  [17-22],  it  is  desirable  to  document  the 
approach  taken  for  this  problem  in  the  present  study. 

The  heat  conduction  analysis  that  relates  the  equivalent 
distributed  source  in  the  form  of  distributed  power  dissipation  to  the 
temperature  distribution  in  the  brain  tissue  has  been  carried  out  by  Guy 
[23,24],  His  approach  was  taken  with  modifications  to  explain  phenomena 
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near  the  interface  of  thermal  discontinuities. 

A.  The  Electromagnetic  Boundary-Value  Problem 

The  basic  approach  has  been  summarized  in  Reference  13  which  is 
included  in  the  Appendix.  In  this  report,  the  discussion  will  be 
restricted  to  the  handling  of  the  numerical  difficulties  in  the  singular 
and  doubly-inf inite  series  Green's  function  employed  in  the  analysis. 

Yaghjian  [17,18]  recently  indicated  that  the  dyadic  Green's 
function  in  the  source  region  was  not  unique  but  rather  dependent  on  the 
geometry  of  the  "principal  volume"  involved.  This  general  statement, 
regardless  of  the  boundary-value  problem  under  consideration,  was 
supported  by  Chen  [25]  for  the  free-space  case.  Their  views  offered  an 
explanation  for  the  apparent  discrepancies  in  the  literature  [19-22]  but 
are  puzzling  on  some  other  aspects.  In  this  report,  no  attempt  is  made 
to  resolve  the  disputes  conclusively  but  merely  to  provide  information 
obtained  in  the  present  analysis  which  could  be  used  by  interested 
researchers. 

The  electric  field  in  the  source  region,  as  shown  in  Figure  2,  can 
be  expressed  in  terms  of  a  dyadic  Green's  function  in  the  following  form 

_E(_r)  =  jyu  lim  /  C(_r,r/)  •  J_(jr')dv’  +  E  (r) ,  (1) 

AV  -*  0  /  V  -  AV  ~C 


where 


£,r’ 

G  (r,r' ) 
G  (r) 
J(r') 

(JO 


positional  vectors  of  the  field  and  source  points, 
the  electrical  dyadic  Green's  function, 
electric  field  at  _r, 
source  current  at  r', 
radian  frequency, 


u  =  permeability  in  space  S, 

V  =  volume  occupied  by  the  volume  current  source, 

AV  -  an  infinitesimal  volume  in  V  enclosing  the 
singularity  of  G^Cjr.j:')  at  £  =  £' ,  and 


E  (r) 
— c  — 


=  a  correction  terra. 


The  space  S  under  consideration  can  be  an  infinite  homogeneous 
space  or  a  homogeneous  space  partially  or  completely  bounded.  The 
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limiting  term  is  often  called  the  principal  value  of  the  integral  since 
its  singular  point  has  been  excluded.  Discussions  on  the  correction 
term  and  the  principal-value  integration  dated  back  to  the  work  of 
Kellogg  [26],  The  work  of  Kellogg  and  Van  Bladel  [27]  implied  that  E 

c 

was  unique  and  independent  of  the  choice  of  AV.  In  treating  the 
principal  integral  of  the  scalar  potential,  Kellogg  explicitly  stated  on 
page  18  that  a  necessary  and  sufficient  condition  for  the  convergence  of 
the  integral  is  that  for  any  small  positive  number  e  there  exists  a 
number  5  >  0  such  that  if  AV  and  AV'  are  any  two  regions  containing 
the  singularity  P  and  contained  in  the  sphere  of  radius  6  about  P, 


(2) 


Kellogg  further  emphasized  that  "the  limit"  is  not  regarded  as 

existing  if  it  is  necessary  to  restrict  the  shape  of  AV  in  order  to 

obtain  a  limit.  The  only  restrictions  on  AV  are  that  it  shall  have  a 

boundary  of  a  certain  degree  of  smoothness  (be  a  regular  region  in  the 

sense  of  Chapter  IV,  Section  8,  p.  100)  that  it  shall  contain  P  in  its 

interior,  and  that  its  maximum  chord  shall  approach  0.  Spheres, 

circular  cylinders,  cubes,  and  pill  boxes  all  have  regular  surfaces,  and 

the  resulting  E  in  Equation  (1)  should  therefore  be  also  independent  of 
— c 

the  shape  of  AV,  even  though  a  vector,  instead  of  a  scalar,  is  involved. 

In  the  present  waveguide  analysis,  AV  is  chosen  to  be  a  rectangular¬ 
sided  cell  whose  dimension  along  the  longitudinal  axis  of  the  waveguide 
approaches  zero  while  dimensions  in  the  other  two  coordinates  remain 
constant.  This  selection  of  principal  volume  AV  appears  to  be  the  same 
as  the  "pill  box"  type  in  Yagjian's  report  [18],  and  is  therefore  in 
agreement  with  all  the  known  views  in  the  literature.  However, 
numerical  experiments  have  not  yet  been  conducted  to  determine  whether 
the  correction  term  E^  is  dependent  upon  the  shape  of  the  principal 
volume  AV. 
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The  self-impedance  matrix  element  Qp  ,  which  corresponds  to  the 

ic| 

limiting  integral  on  the  right-hand  side  of  Equation  (1),  should  hold 

the  answer  to  the  present  puzzle.  Since  the  electric  field  E_(jr)  is 

unique,  the  two  terms  on  the  right-hand  side  of  Equation  (1)  must  add  up 

to  yield  the  correct  value  for  EOr).  One  can  therefore  construct 

various  principal  volumes  to  see  whether  a  dependence  on  the  principal 

volume  exists  for  the  correction  term  E  ,  and,  if  it  does,  whether  a 

— c 

generalization  as  claimed  by  Yagjian  [18]  is  valid.  Even  though  no 
formal  investigation  has  yet  been  made  in  this  respect,  an  examination 
of  the  mathematical  expressions  in  the  waveguide  case  [13],  the  free 
space  case  [28],  and  the  strict  convergence  criterion  of  Equation  (2) 
indicates  that  the  issue  is  far  from  being  resolved. 

The  self-impedance  element  Q^,  for  p  “  1,  can  be  evaluated  by  a 
partial  summation  technique  as  follows: 

Q11  =  T^  +  s11 

%  £q  *£q  ’  (3) 


where 

„11 


11  «  _J_  'O  V  /nit  \  2  ^  .2  4ab 

eq  2abk=  Zw  2^  b  )  nm  nirmir 
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and  all  the  notations  are  the  same  as  defined  in  Reference  13. 

The  term  in  Equation  (3)  converges  rapidly  but  the  term 
converges  extremely  slow.  Using  the  summation  formulas  in  [13]  one 
obtains 
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and,  for  2  -  k2  >  0  , 
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Similarly,  expressions  for  Q  can  be  obtained  by  exchanging 
between  m  and  n,  between  AX^  and  AY^,  between  X^  and  Y^,  between  X^  and 
Y  ,  and  between  a  and  b  in  Equations  (3)  through  (8).  The  expressions 
for  Q  are 

'  iq 


+  S 


33 

iq 


(9) 


11 


where 


33 

iq 


2abk2 


ro 


m 


m=0 


4a  b 
nimn 


miAX. 


(sin 


2a 


mi  AY 


sin 


~2T>  jir  coslk„DD-D)) 


nm 


nm 


eXp(jknm  AZi/2)  sin 


nTrX  nirX^  mil  y 

sin  - -  sin  — — 

a  a  b 


irn  Y 


sin 


CO) 


T33  =  T33  +  T33 

D  a  *q  b  iq 


CD 


In  Equation  (11), 
_  1 


AY, 


aiq  =  4k^  j  b(Yq  +  Y9,  +  -T')  ~  ^Yq  +  Yf_ 


D, 

9  ; 


r  AV  i 

Id  AXj>  1  AX 

•|a(Xq  +  \  +  1“)  ~  J<Xq  +  -  -±) 


+  V  [1 


D 

a  q 


AX. 


(Xa  “  Xf  +  1  ~  n*  ( 1  -  r’  (X  -  X  -  i)  ] 


a  q 


)i(,  (1 


where 


AY? 

f°r  Yq  “  Yn  +  —  D 
otherwise  » 

fot  Y,  -  Y<  -  -r  * 0 

otherwise 

AX, 

for  X  -  X,  -  !  0 

q  >-  2 

otherwise 


c  AX 

for  Xq  -  X, +  2  0 

otherwise 


(13) 


12 


Also  in  Equation  (11), 
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It  appears  that  numerical  experiments  can  be  fruitfully  conducted 

by  letting  AX^  and  AY „  vanish  as  rapidly  as  AZ^  and  then  observe  what 

correction  term  E  (r)  should  be  used.  Short  of  a  laboratory  experiment, 

— c  — 

these  tests  should  be  very  useful  in  resolving  the  existing 
discrepancies  in  the  literature. 

It  is  also  interesting  to  note  that  a  similar  singular  property  has 
also  been  observed  in  Green's  function  representation  of  an  infinite 
planar  array  with  wire  antenna  elements.  It  was  observed  that  the 
Green's  function  has  a  singular  point  associated  with  the  component  of 
the  equivalent  source  current  in  the  direction  perpendicular  to  the 
plane  of  the  array  [29,30].  The  computer  program  yielded  good  results 
when  the  array  was  made  of  horizontal  dipoles  (dipoles  parallel  to  the 
array  plane),  but' it  "blew  up"  when  dealing  with  vertical  dipoles 
(dipoles  perpendicular  to  the  array  plane).  Since  the  planar  array 
problem  is  intimately  related  to  the  rectangular  waveguide  problem,  this 
phenomenon  seems  to  support  the  result 


in  References  20,  21,  22  and  dispute  Yagjlan's  [18]  result  which  states 
that 


E  -  d/3, 

— c 

where  a  cube  is  chosen  as  the  principal  volume. 


(16) 
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B.  The  Temperature  Distribution  Inside  a  Biological  Body  Under 
Microwave  Excitation 

.  The  density  of  power  dissipated  in  a  biological  body  is 

pd(r)  =  UE-(r)  |E(r)|,  (17) 

where  e"  is  the  imaginary  part  of  the  complex  permittivity.  The  density 
of  power  dissipated,  pd  in  Equation  (17),  can  be  considered  as  a 
distributed  heat  source  inside  the  biological  body.  The  temperature 
distribution  as  a  function  of  time,  space, and  initial  conditions  can 
then  be  computed  from  purely  thermal  considerations.  The  problem  can  be 
formulated  by  the  heat  equation  [31]  as  follows 

3T  2  2  2 

—  =  a  V  T  +  4-rra  S  .  (id) 

where  T  ■  temperature  distribution  function, 
t  *  time, 

V2  *  Laplacian  operation, 

a  =  K/pc,  K  being  the  heat  conductivity,  p  the  density,  and 
c  the  specific  heat  of  the  medium, 

S  -  P  /4  nK ,  and 
d 

P,  =  heat  source  distribution, 
a 

Both  T  and  P,  in  Equation  (18)  are  functions  of  time  t  and  spatial 
d 

coordinate  r.  P  represents  the  distributed  source  which  is  a  known 
quantity  computed  from  Equation  (17).  The  solution  of  Equation  (18) 
together  with  all  the  initial  and  boundary  conditions  has  been  obtained 
for  several  elementary  geometries  [31,32].  For  complex  geometries, 
numerical  matrix  method  such  as  the  finite-difference  method  can  be 
used.  For  simple  cases,  an  approximate  formula  used  by  Guy  [23]  is  very 
convenient.  The  short-term  temperature  rise  AT  is  approximately 
proportional  to  the  power  dissipation  at  the  point  of  interest,  namely, 
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p 


d 


0  CAT, 


(19) 


where  p  and  c  are  as  defined  in  Equation  (18)  and  a  is  a  coefficient 
related  to  the  exposure  time.  Based  on  Equation  (19),  the  relative 
short-time  temperature  distribution  is  approximately  the  same  as  the 
density  of  power  dissipation.  Guy  did  not  indicate  the  origin  nor  the 
rationale  for  Equation  (19).  However,  Equation  (19)  can  be  considered 

to  be  a  special  case  of  Equation  (18)  at  a  short  period  after  t  =  0. 

2  2 
Since  7  T  *  0  at  t  =*  0-,  as  was  the  usual  case,  we  have  V  T  =  0  at  T  = 

0+,  and  Equation  (18)  can  be  approximated  by  Equation  (19). 

Equation  (19)  is  often  satisfactory  when  applied  to  the  interior 
region  of  a  homogeneous  portion  of  a  dielectric  body  of  low  thermal 
-'(inductivity.  It  fails  to  provide  useful  results  near  the  edge  of  the 
dielectric  body  or  at  the  interface  of  two  different  types  of 
dielectrics,  as  was  noted  by  Guy.  As  shown  in  Figure  3,  Guy's 
measurement  of  temperature  prof!  of  a  block,  formed  with  muscle  in 
Region  1  and  fat  in  Region  2,  exhibited  gross  deviation  at  the  interface 
f rora  his  -alculated  power  dissipation  curve.  Thus,  Equation  (19) 
appeared  to  fail  in  the  neighborhood  of  the  interface  between  the  fat 
and  the  muscle  tissues.  Guy  then  proceeded  to  correct  this  discrepancy 
b v  enforcing  the  continuity  of  temperature  at  the  interface.  Although 
his  approach  appears  reasonable  on  a  qualitative  basis,  many  questions 
remain  to  be  answered. 

A  more  rigorous  appraoch  was  taken  in  the  present  study  by  using 
the  diffusion  Equation  (18).  We  first  investigated  the  fat-muscle 
block  of  Figure  A,  and  was  able  to  obtain  good  numerical  results. 
However,  we  have  found  that  the  temperature  distribution  is  quite  time 
dependent,  being  fastly  varying  in  the  first  five  seconds.  Thus,  the 
temperature  profile  can  not  be  simply  established  by  recording  within, 
say,  five  seconds  after  the  heating. 

We  first  examine  the  case  in  which  the  muscle  is  held  at  a  uniform 
temperature  Tq  and  the  fat  is  held  at  0°k  at  time  t  =  0-  as  shown  in 
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Figure  4.  In  the  absence  of  a  heating  source,  the  temperature 
distribution  is  [32] 
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and  k^,  C^,  and  k^,  are  the  heat  conductivity  and  specific  heat  for 

Regions  1  and  2,  respectively. 

The  temperature  profile  along  the  x  axis  changes  rapidly  during  the 
first  five  seconds,  as  shown  in  Figure  5.  As  a  result,  merely  recording 
the  temperature  as  soon  as  possible  (for  example,  Guy  did  it  in  less 
than  five  seconds)  does  not  guarantee  consistent  and  meaningful  results. 
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mpernture  profile  along  a  composite  block  of  fat  and  muscle. 


Equation  (20)  can  also  be  employed  to  compute  the  approximate 
temperature  for  the  case  in  which  the  initial  temperature  distribution 
is  not  uniform.  An  example  is  shown  in  Figure  3  for  the  fat-muscle 
block  studied  by  Guy.  The  relative  temperature  distribution  can  be 
approximated  by  modifying  the  temperature  expression  of  Equation  (20) 
with  the  nonuniform  dissipated  power  distribution  P,.  As  can  be  seen, 
the  temperature  peak  near  the  interface  on  the  muscle  side  shifts  along 
the  x  axis  during  the  first  10  seconds.  This  observation  is  important 
in  understanding  the  thermographic  data  obtained  in  the  present 
research,  in  which  an  apparent  hot  spot  near  the  edge  of  a  tissue  has 
also  been  observed. 

It  is  difficult  to  solve  Equation  (18)  for  a  biological  tissue  of 
more  complex  geometry  other  than  slabs,  spheres,  ellipsoids,  and 
cylinders.  However,  it  has  been  observed  that  the  simple  approximation 
of  Equation  (19)  is  valid  in  regions  far  from  sharp  temperature 
gradients  such  as  the  edge  of  or  interface  in  a  tissue. 
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SECTION  III 


COMPARISON  WITH  THERMOGRAPHIC  MEASUREMENTS 

Two  types  of  thermographic  measurements  were  made:  one  using 
thermographic,  paper  and  the  other  using  thermographic  camera.  The 
thermograhic  paper  measures  the  peak  temperature  during  and  after  the 
exposure  time  but  does  not  record  accurate  data.  Furthermore,  it  is 
always  difficult  to  place  a  thermographic  paper  at  the  proper  location 
without  affecting  the  experiment.  The  thermographic  camera  is  more 
accurate  hut  is  also  difficult  to  look  at  the  inside  of  a  tissue.  As 
has  been  discussed  in  the  preceding  section,  the  temperature 
distribution  varies  with  time.  Merely  taking  a  quick  picture 
immediately  after  heat  exposure  does  not  always  yield  consistent  data. 
Consequently,  great  care  must  be  taken  when  comparing  the  calculated  and 
measured  thermal  data.  Measurements  were  made  at  WRAIR,  even  though 
some  preliminary  tests  were  performed  at  Georgia  Tech.  All  the  cases 
studied  were  rectangular-sided  blocks  which  are  convenient  for  the  book¬ 
keeping  of  the  location  of  the  field  points. 

Figure  6  shows  a  comparison  between  the  calculated  dissipated  power 
intensity  and  the  temperature  profile  recorded  with  thermographic  paper 
for  Case  C  in  Reference  13.  Although  a  general  agreement  can  be 
observed  in  the  comparison,  the  thermographic  measurement  was  difficult 
to  calibrate  and  repeat.  The  dielectric  model  used  had  a  high  water 
content  which  generated  hot  steam  during  the  later  stage  of  the 
exposure.  A  power  level  of  15  to  50  watts  was  used,  and  exposure  times 
of  15  to  25  minutes  were  tried.  Usually  four  or  more  recordings  of 
different  exposure  times  must  be  tried  before  a  recording  of  good 
quality  can  be  obtained. 

Figure  7  shows  the  physical  configuration  of  a  dielectric  (Stycast 
Hi-KIb)  block  inside  a  WR284  waveguide.  The  block  is  divided  into  54 
identical  cells  which  were  numbered  as  shown  in  the  figure.  The 
calculated  dissipated  power  intensity  at  the  center  of  each  cell  on  the 
top  layer  is  plotted  in  Figure  8  in  a  perceptive  view.  Two 
temperature  profiles  measured  with  a  heat-sensitive  camera  are  shown  in 
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Figure  6.  Calculated  distribution  of  dissipated  power  in  comparison  with  thermographic  mea 


A  dielectric  block  inside  a  rectangular  WR284  waveguide. 


Figure  8.  Di.st 
pass 
bloc 


Figure  9  for  comparison.  Satisfactory  agreement  is  observed  except  for 
a  discrepancy  on  the  far  end  of  the  block  where  the  calculated  power 
dissipation  continuously  rises  toward  the  edge  of  the  block  while  the 
measured  temperature  shows  a  drop  near  the  edge.  This  disagreement  can 
be  explained  in  the  following  manner. 

As  has  been  discussed  in  the  preceding  section,  the  temperature 
distribution  is  approximately  proportional  to  the  dissipated  power 
intensity  only  for  a  brief  period  after  the  exposure  and  only  for 
regions  away  from  the  surfaces  with  discontinuity  of  thermal 
coefficients.  The  drop-off  of  the  temperature  near  the  distal  (load) 
end  of  the  block  can  be  explained  in  the  same  manner  it  was  in  Figures  3 
and  5.  The  calculated  power  dissipation  corresponding  to  (a)  of  Figure 
9  is  plotted  in  Figure  10  in  the  same  fashion.  Good  numerical 
convergence  is  indicated  in  computations  using  24,  54,  and  63  cells. 

The  63-cell  calculation  employed  a  configuration  of  unequal  cell  volumes 
to  obtain  great  details  in  the  questioned  distal  end  near  z  =  2  cm.  As 
can  be  seen,  there  is  a  definite  rise  in  power  dissipation  toward  the 
edge  of  the  block.  Temperature  distribution  is  calculated  using  thermal 
conductivities  of  0.0084  and  0.0008  watts/cm-°K  for  specific  heats  of 

3 

4.L86  and  0.1  joules/cm  -°K  for  the  dielectric,  and  the  air, 
respectively.  The  resulting  temperature  fall-off,  however,  was  much 
closer  to  the  edge  than  the  measured  data  in  Figure  9(a).  We  have 
found,  however,  that  the  temperature  peak  can  be  shifted  away  by  as  much 
as  0.4  cm  from  the  edge  by  using  different  thermal  coefficients  in  the 
computation.  Therefore,  the  calculation  could  be  improved  by  using  more 
reliable  thermal  coefficients  and  by  a  more  refined  thermal  modeling 
approach. 


(a)  Along  the  waveguide  axis 


(b)  Perpendicular  to  the  waveguide  axis 


Figure  9.  Thermographic  measurements  on  the  top  surface  of  the 
block  sample  of  Figure  7. 


SECTION  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

An  analytical  tool  has  been  developed  to  compute  the 
e lectromagnetic  fields  inside  a  three-dimensional  arbitrarily-shaped 
heterogeneous  dielectric  or  biological  body  inside  a  rectangular 
waveguide.  The  profile  of  dissipated  power  intensity,  and  hence  the 
temperature,  can  be  determined  with  satisfactory  accuracy.  Immediate 
application  of  this  technique  will  be  to  provide  guidance  for  the  design 
of  microwave  applicators  for  rapid  inactivation  of  enzymes.  The 
techniques  involved  in  the  present  analysis  represent  advances  in 
microwave  and  electromagnetic  theory.  Their  impact  could  be  extended 
into  other  physical  and  mathematical  sciences. 

Numerical  tests  have  been  conducted  for  rectangular-sided 
homogeneous  bodies  only,  even  though  the  computer  program  was  written 
for  a  general  three-dimensional  arbitrarily-shaped  heterogeneous  body. 
The  computer  program  was  checked  by  numerical  convergence  tests  and 
conservation  of  energy,  as  well  as  by  direct  comparison  with  measured 
reflection  and  thermographic  data.  For  dielectric  bodies  of  low  to 
medium  dielectric  constants  (less  than,  say,  lb),  excellent  results 
were  obtained.  The  quality  of  the  results  degraded  gradually  with 
increasing  dielectric  constants.  However,  even  for  high  dielectric 
constants  of  about  50,  fairly  good  results  were  yielded  by  the  computer 
pro  -ram. 

i’urther  research  should  include  test  cases  of  biological  bodies  of 
greater  complexity,  such  as  a  two  layer  ellipsoid  or  even  a  full  rat 
configuration.  In  addition,  the  waveguide  structure  can  be  modified  to 
represent  a  more  practical  applicator  such  as  a  waveguide  with  a 
shorting  plate  at  the  end.  Ttiese  tasks  have  been  Initiated  during  the 
present  grant  period  but  could  not  be  completed  with  the  limited  funds 
ava 1 lable. 

It  Is  also  recommended  that  a  more  accurate  theoretical  approach  be 
established  to  relate  the  distributed  power  dissipation  and  the 
temperature  profile  Inside  the  biological  body.  Extensive  numerical 
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data  should  be  generated  to  guide  the  applicator  design.  However, 
improvements  of  the  existing  algorithm  must  be  made  before  data 
generation,  since  the  computer  runs  are  quite  costly. 
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Analysis  of  a  Three-Dimensional  Arbitrarily 
Shaped  Dielectric  or  Biological  Body 
Inside  a  Rectangular  Waveguide 

JOHNSON  J.  H.  WANG,  member,  ieee 


Abstract — This  paper  presents  a  method  for  the  analysis  of  three-dimen¬ 
sional  arbitrarily  shaped  dielectric  obstacles  inside  a  rectangular  ware- 
guide.  The  numerical  computation  involves  a  dyadic  Green’s  function 
containing  a  double  infinite  series,  which  is  evaluated  by  a  partial  summa¬ 
tion  technique. 

I.  Introduction 

WAVEGUIDE  obstacles  and  discontinuities,  includ¬ 
ing  the  dielectric  type  to  be  discussed  in  this  paper, 
are  long  standing  problems  in  electromagnetic  theory. 
Many  of  them,  essentially  two  dimensional,  have  been 
solved  and  were  summarized  by  Marcuvitz  [1],  The  gen¬ 
eral  three-dimensional  discontinuity  problems,  however, 
remain  unsolved  in  spite  of  the  advent  of  modem  high¬ 
speed  digital  computers  and  the  method  of  moments  [2] 
which  permit  treatment  of  problems  not  solvable  by  exact 
methods.  Upon  reviewing  the  status  of  numerical  tech¬ 
niques  for  passive  microwave  devices.  Silvester  and 
Csendes  [3]  observed  in  1974,  “not  a  single  truly  three-di¬ 
mensional  solution  has  been  published”  for  waveguide 
discontinuity  problems.  This  statement  is  apparently  still 
valid  today. 

This  lack  of  published  research  activities  in  three-dimen¬ 
sional  waveguide  discontinuities  has  been  in  many  cir¬ 
cumstances  due  to  the  deficiencies  of  the  Green's  func¬ 
tions  in  the  waveguide  region.  A  dyadic  Green's  function 
for  the  rectangular  waveguide  was  presented  by  Tai  in 
1972  [4]  and  later  revised  by  the  same  author  in  1973  [5]. 
Tai's  expression  includes  a  double  infinite  series  summing 
over  the  contributions  from  all  the  individual  waveguide 
modes.  Recalling  the  simple  expression  of  the  free  space 
Green’s  function,  one  immediately  recognizes  the  greater 
complexity  in  the  waveguide  case. 

This  paper  presents  a  successful  use  of  the  dyadic 
Green's  function  in  the  analysis  of  the  electromagnetic 
problem  of  a  three-dimensional  arbitrarily  shaped  dielec¬ 
tric  or  biological  body  inside  a  rectangular  waveguide.  In 
the  process,  the  extremely  slow  convergence  of  the  double 
infinite  series  in  the  Green's  function  had  to  be  modified 
by  means  of  a  partial  summation  technique.  The  im¬ 
mediate  application  of  this  new  analytical  technique  is  in 
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Fig.  1.  A  three-dimensional  arbitrarily  shaped  heterogeneous  dielectric 
body  illuminated  inside  a  rectangular  waveguide 

microwave  waveguide  enzyme  inactivation  in  neuro-chem¬ 
ical  research  [6].  Extension  of  this  technique  to  highly 
conductive  obstacles  is  feasible  but  may  require  modifica¬ 
tion  of  the  volume  integral  into  a  surface  expression. 

II.  The  Integral  Equation  and  the  Dyadic 
Green's  Function 

The  problem  to  be  considered  is  shown  in  Fig.  1.  in 
which  a  three-dimensional  arbitrarily  shaped  heteroge¬ 
neous  dielectric  or  biological  body  is  electromagnetically 
illuminated  in  a  rectangular  waveguide.  The  dielectric 
body  has  a  permittivity  distribution  of  e(r),  where  r  is  a 
positional  vector.  Outside  the  volume  V  occupied  by  the 
dielectric  body,  the  perrmttivity  is  homogeneous  and  is 
denoted  by  Free  space  permeability  ^  is  assumed  for 
both  the  dielectric  body  and  the  medium  outside  U. 

The  time  function  where  i  and  m  are  time  and 

radian  frequency,  is  used  in  all  the  equations  in  this  paper 
for  ready  comparison  with  Tai's  work  [4],  [5].  Since  the 
e'*"  convention  is  perhaps  more  widely  used,  a  comment 
on  the  conversion  of  the  equations  to  this  convention  is 
justified  at  this  point.  For  the  e/"'  convention,  one  merely 
changes  to  —j  all  the / s  appearing  in  the  equations  in  this 
paper. 

In  Fig.  1,  £(r)  denotes  the  electric  field  intensity  at  r 
and  E'(r)  denotes  the  field  intensity  at  r  with  the  dielec¬ 
tric  body  replaced  by  the  medium  The  scattered  field  is 
defined  as 


£T(r)  =  E(r)-E'(r).  (1) 

The  volume  equivalence  principle  can  be  shown  to  be 
valid  in  the  bounded  as  well  as  the  unbounded  space.  As  a 
result,  the  problem  as  depicted  in  Fig.  I  is  equivalent 
everywhere  to  a  homogeneous  waveguide  with  ^  and 
with  volume  current  density 

/(/■)=  -ywfefr)-*,  j£(r).  (2) 
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It  follows  from  the  theorem  of  superposition  that  the 
electric  field  radiated  from  J{r)  is  equal  to  E’(r).  There¬ 
fore, 


E’(r)=j^f  GAryyJWdV  (3) 

where  Ge  is  the  dyadic  Green’s  function  of  the  electric 
type  and 

G>.'0=^(r,rO--!riiS(r,r')  (4) 

where  S  denotes  a  three-dimensional  Dirac  delta  function, 
z  denotes  a  unit  vector  along  c,  kK  =(jVe,  ft 0  ,  and  is 
defined  below.  The  term  involving  S  had  not  been  in¬ 
cluded  in  Tai’s  earlier  work  [4]  until  1973  (5]  and  has  been 
a  subject  of  recent  discussions  [5],  [7],  [8].  was  not 
explicitly  expressed  in  the  literature  except  by  Rahmat- 
Samii  [7],  However,  Rahmat-Sarrui’s  expression  has  a 
number  of  errors  in  the  print.  For  the  clanty  of  the 
present  discussion  and  the  convenience  of  future  refer¬ 
ences,  it  is  desirable  to  present  it  in  the  following  long 
form: 


Ge0(r.r')  = 


/  XI  X  f  f 

J  y  y  iQn  Qmc/^u-^i 
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Equations  ( 1 )— (4)  can  be  manipulated  to  yield  the 
following  integral  equation: 

c  -  J(r) 

/uififl  /  G^(r.r')J(r)  d\  '+ - r - - 

Jv  ywi  Ctrl  -  1 1 

J.  ( r): 

- - =  -  E  in  (61 
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where  J t  is  the  c  component  of  J  The  unknown  J  in  (6i 
can  then  be  solved  by  the  method  of  moments  [2]  which 
transforms  the  integral  equation  into  a  set  of  linear  equa¬ 
tions  readily  solvable  by  means  of  a  digital  computer. 
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III.  SOLUTION  BY  THE  METHOD  Ot  MOMENTS 

Although  there  exists  a  number  of  methods  by  which  an 
integral  equation  can  be  solved  numerically,  the  complex¬ 
ity  of  three-dimensional  geometry  can  hardly  tolerate  fur¬ 
ther  complication  in  the  computational  process.  Even  in 
the  much  simpler  case  of  plane  wave  incidence  and  un¬ 
bounded  free  space,  only  point  matching  together  with 
rectangularly  sided  cells  has  been  attempted  for  the 
volume  type  of  integral  equation  [9],  Fortunately,  this 
unsophisticated  process  has  been  found  to  be  capable  of 
producing  good  numerical  results.  Thus  point  matching 
with  rectangularly  sided  cells  is  employed  in  the  present 
analysis. 

The  volume  V  occupied  by  the  dielectric  body  is  first 
divided  into  L  equal  rectangular-sided  cells  AT, -AT,, 
each  of  which  has  constant  dimensions  Ax.  Ay.  and  Ac. 
The  electric  field,  assumed  to  be  uniform  inside  the  /th 
cell,  is  designated  as  £(/■,),  where  r,  represents  the  center 
of  the  /th  cell.  The  equivalent  current  in  (2)  can  then  be 
expressed  as 

L  1 

y(r)- 2  2  Win  (?) 

/-I  k -  I 

where  B*  is  a  basis  function  defined  as 

£,*(/■)  k=>  1.2.3  or  v.v.c  (8) 

uk  in  (8)  denotes  a  unit  vector,  and 
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The  weighting  function  is  defined  as 

Wr{r)  =  8{r-rq)up,  p- 1,2,3.  (10) 

The  scalar  product  between  /  and  g  is  defined  as 

(/*)=  ffgdv.  (11) 

Jv 

We  can  generate  a  set  of  linear  equations  by  first 

substituting  (7)  into  (6)  and  then  performing  a  scalar 

product  on  the  resulting  equation  with  the  weighting  func¬ 
tion  of  (10)  for  p=  1,2,3  and  q=*  1,-  •  ■  ,L.  This  moment 
generating  procedure  leads  to  the  following  set  of  equa¬ 
tions: 

2  2  J,kAtf=  Cpq,  p=  1,2,3;  q=\,-  -  ■  ,L  (12) 

*-i  /-i 

where 
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S£  in  ( 14)  is  the  ICronecker  delta,  being  1  when  p  =  q  and  0 
otherwise,  and 

Qfqk  =  f  G$(rrr')dV  (15) 

-'•sr, 

where  G%  is  the  (p.k)  component  of  the  dyad  G. o  The 
integration  in  (15)  can  be  earned  out  to  yield 
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and  /  in  (16)  depends  on  the  index  pk.  For pk  being  (1,3), 
(2, 3),  (3, 1 ).  (3,2),  where  ( 1, 2, 3)  corresponds  to  (.x.y, j),  we 
have 


For  other  pk  indices,  we  have 
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Fig.  2.  Comparison  of  convergence  between  direct  truncation  and  par¬ 
tial  summation. 


As  can  be  seen  in  ( 14)— (20).  the  matnx  element  Aff 
involves  a  double  infinite  senes  which  does  not  converge 
rapidly  unless  |2?-2,|  is  large.  When  p  =  q  as  is  the  case 
for  the  self  cells  which  are  the  diagonal  matrix  elements 
|z?  — 2,|  =  0  and  the  series  convergence  is  extremely  slow. 
As  an  example.  Fig.  2  shows  that  no  sign  of  convergence 
is  exhibited  even  after  140x  140  terms  are  used  in  the 
summation.  Thus  the  computation  of  the  matnees  can  not 
be  carried  out  by  a  simple  truncation  of  the  series. 

This  computational  difficulty  can  be  surmounted  by 
direct  summation  of  the  part  of  the  infinite  series  which 
contains  terms  nonvanishing  with  increasing  .V/  and  A. 
This  partial  summation  transforms  a  matnx  element  of  a 
double  infinite  series  into  a  single  infinite  senes  or  even 
into  an  expression  of  closed  form,  which  can  be  truncated 
for  computation.  Fig.  2  shows  that  convergence  is 
achieved  with  about  20x20  terms  for  the  partial-summa¬ 
tion  technique,  which  is  significant  improvement  over  the 
direct  truncation  method.  For  most  off-diagonal  matrix 
elements,  direct  truncation  is  satisfactory  because  of  the 
predominant  influence  of  the  exponential  terms  with  the 
argument  of  jkmn\zq  -  z,\  or  jkmnl:t/2  as  shown  in  (19) 
and  (20).  As  a  result  of  these  exponential  terms  which 
rapidly  decrease  with  m  and  n,  off-diagonal  matrices  with 
nonvanishing  |a,  —  j  |  can  be  computed  with  a  finite  senes 
truncated  according  to  a  precision  entenon  established  by 
the  value  of  jkmn\zq  -  2,|.  Depending  on  the  value  of  j;^- 
2, j,  approximately  16x8  up  to  21  x  12  terms  were  used  for 
m  and  n  in  the  examples  reported  in  this  paper. 

The  partial  summation  technique  is  tedious  but 
straightforward.  The  portion  of  the  infinite  senes  which 
consists  of  terms  slowly  convergent  with  m  and  n  are 
summed  by  means  of  the  following  formulas  [10): 
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Itrnr  <  x  <(2m  +2)v,  a  being  noninteger. 

The  explicit  expressions  for  the  matrix  elements  after 
partial  summation  are  too  complex  to  be  presented  in  this 
paper  and  the  interested  readers  are  referred  to  {11}  for 
further  details.  Their  computation  on  a  digital  computer  is 
efficient  since  only  logical  expressions  and  elementary 
algebraic  and  transcendental  functions  are  involved. 

IV.  Numerical  Examples  and  Supporting 
Measurements 

A  method  for  analyzing  the  electromagnetic  problem  of 
a  three-dimensional  arbitrarily  shaped  body  inside  a  rect¬ 
angular  waveguide  has  been  presented.  A  Fortran  IV 
computer  program  was  written  for  the  numerical  testing 
of  this  approach.  Measurements  were  also  conducted  to 
compare  with  the  theoretical  prediction.  Several  cases 
were  studied  and  satisfactory  results  have  been  obtained. 
Validation  of  the  theory  and  the  computer  program  was 
achieved  with  numerical  convergence  tests  as  well  as 
direct  comparison  with  measured  data,  which  included 
transmission-reflection  characteristics  and  thermographic 
heating  patterns. 

Three  cases,  as  shown  in  Fig.  3,  are  presented  in  this 
paper.  All  of  the  test  cases  consist  of  homogeneous  dielec¬ 
tric  bodies  with  rectangular  side  walls  aligned  with  the 
waveguide  walls.  This  choice  of  geometry  conformal  to 
the  Cartesian  coordinates  is  mainly  for  the  sake  of  sim¬ 
plicity  in  data  management  and  should  not  result  in  any 
significant  loss  of  generality  as  was  noted  in  the  free  space 
case  (9).  For  highly  conductive  obstacles,  the  surface 
curvature  of  the  obstacle  plays  a  more  important  role  and 
therefore  should  be  treated  with  greater  discretion. 

It  was  noted  that  the  linear  cell  dimensions  should  be 
A/2  (A  being  the  free  space  wavelength  divided  by  the 
square  root  of  the  dielectric  constant)  or  less  in  order  to 
yield  accurate  data.  This  observation  had  been  reported  in 
the  free  space  case  studied  by  Livesay  and  Chen  [9]  and 
later  by  Hagmann  el  al.  [12],  Fig.  4  shows  good  agree¬ 
ment  in  the  reflection  and  transmission  properties  of  Case 
B  between  a  12-cell  calculation  and  the  measured  data 
using  a  model  made  of  silica  compound.  A  12-cell  calcula¬ 
tion  from  Case  A.  being  also  a  case  of  low  dielectric 
constant,  yields  a  power  reflection  coefficient  of  0.114  at 
2.65  GHz.  dropping  down  to  0.035  at  3.5  GHz,  which  was 
verified  experimentally  with  a  paraffin  wax  model. 

While  12  cells  are  sufficient  for  the  calculation  of  Cases 
A  and  B,  many  more  cells  are  needed  for  Case  C,  which 
has  a  high  dielectric  constant.  Case  C  was  originally 
intended  to  be  a  phantom  model  for  simulating  muscle 
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Fig.  3.  Configurations  of  the  three  cases  studied. 
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Fig.  4.  Comparison  between  calculated  and  measured  reflection  and 
transmission  charactensucs  for  Case  B  of  Fig.  3  which  is  of  low 
dielectric  constant. 
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Fig.  5.  Comparison  between  measurements  and  calculations  of  various 
nuhibers  of  cells  for  Case  C  of  a  high  dielectric  constant 

tissue.  The  model  was  made  by  mixing  water,  powdered 
polyethylene,  and  "super  stuff.”  a  modeling  compound 
manufactured  by  Wham-0  Co.,  San  Gabriel.  CA.  The 
dielectric  constant  and  loss  tangent  were  then  measured  ai 
various  frequencies.  There  were  difficulties  in  achieving 
and  maintaining  the  desired  dielectric  properties  of  the 
phantom  model  and.  as  a  result,  only  the  repeatable 
measured  data  are  shown  in  Fig.  5.  in  which  comparisons 
are  also  made  for  three  calculations  using  different  num¬ 
bers  of  subvolume  cells.  The  measured  data  in  Fig.  5 
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NUMBER  OF  CELLS 

Fig.  6.  Convergence  of  transmission  and  reflection  properties  for 
C  with  high  and  medium  dielectric  constants  at  2.8  GHz. 


try  of  Case  C  with  two  values  of  complex  dielectric 
constant.  In  Fig.  7,  the  field  distribution  inside  the  dielec¬ 
tric  body  is  compared  for  various  numbers  of  subvolume 
cells  used  in  the  calculation.  It  has  been  noted  that  the 
convergence  of  field  intensity  is  related  to  the  profile  of 
the  field.  When  the  field  vanes  slowly  with  the  coordi¬ 
nates,  convergence  is  rapid.  Fig.  7  presents  typical  situa¬ 
tions  with  moderately  varying  fields.  The  favorable  in¬ 
fluence  of  the  lower  dielectric  constant  on  the  conver¬ 
gence  of  the  field  distribution  is  clearly  demonstrated. 

Dielectric  bodies  of  high  dielectric  constant  require  not 
only  a  larger  number  of  cells  but  also  a  greater  number  of 
terms  in  the  Green’s  function  series.  This  is  due  to  the  fact 
that  the  distance  between  adjacent  cells  is  small  because 
of  the  smaller  size  of  the  cells.  The  attenuation  of  higher 
order  modes,  being  independent  of  the  dielectric  constant, 
decreases  when  the  distances  between  adjacent  cells 
centers  are  shortened.  As  a  result,  the  calculation  cost  for 
cases  of  high  dielectric  constants  increases  rapidly  with  an 
increase  in  dielectric  constant. 
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Fig.  7.  Convergence  of  field  distribution  for  Case  C  at  2.8  GHz  with 
high  and  medium  dielectric  constants. 

agree  reasonably  with  the  calculation.  The  12-cell  config¬ 
uration  is  obviously  overly  crude,  yet  the  results  are  in 
gross  agreement  with  those  of  finer  configurations.  Note 
that  in  the  45-cell  calculation  the  >>  dimension  of  each  cell 
is  0.8833  cm,  which  is  about  0.52A  at  2.5  GHz  and  0.68A  at 
3.1  GHz.  C  nvergence  becomes  more  rapid  when  the 
linear  dimensions  of  the  subvolume  cells  decrease  to  A/2 
or  less. 

The  rapidity  of  convergence  of  the  present  model  ap¬ 
proach  is  further  illustrated  in  Figs.  6  and  7.  Fig.  6  shows 
the  calculated  power  reflection  and  transmission 
coefficients  versus  the  number  of  cells  used  for  the  geome- 


V.  Concluding  Remarks 

A  general  three-dimensional  waveguide  dielectric  ob¬ 
stacle  has  been  successfully  treated  by  employing  the 
moment  method  on  an  integral  equation  involving  a  dy¬ 
adic  Green’s  function.  This  general  method  can  be  ap¬ 
plied  to  a  number  of  waveguide  problems.  An  immediate 
extension  of  this  technique  would  be  to  ferromagnetic 
obstacles.  For  highly  conductive  obstacles  a  surface-type 
Green’s  function  may  be  more  desirable. 
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